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ABSTRACT 


Two Extended Kalman Filter algorithms for air-to-air 
passive ranging are proposed, and examined by computer 
simulation. One algorithm uses only bearing observations 
while the other uses both bearing and elevation angles. 

Both are tested using a flat-Earth model and also using a 
spherical-Earth model where the benefit of a simple correc- 
tion for the curvature-of-the-Earth effect on elevation angle 
is examined. The effects of varied angle measurement 
precision and of varied range estimate accuracy are inves- 
tigated. Both filters are found to perform acceptably under 
most conditions, with the bearing-elevation algorithm 


generally providing superior results. 
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The problem of Fleet Air Defense and air defense in 
general is one covering increasingly large geographical 
areas due to the greater speeds, ranges, and missile launch 
envelopes of potentially hostile airborne attack platforms. 
Concurrent with the expansion of the physical size of the 
problem has been an increase in its complexity particularly 
in the area of Electronic Warfare (EW). 

A common effect of EW is the denial of range and range 
rate information by electronic jamming of the observer's 
radar system. Since target range and range rate must 
generally be known or at least estimated for a fighter to 
successfully employ its missile system, the fighter must 
have some alternate means of obtaining that information. 

An approach to this problem is called Passive Ranging 
(PR) or pasSive Target Motion Analysis (TMA). This might 
involve one or more observers using a "receive-only" radar 
mode or perhaps a highly directional infrared (IR) sensor. 
It could entail such things as: (a) triangulation of 
bearing observations from more than one observer, (b) use of 
Ome, Or more known target parameters (with continued bearing 
observations) to compute the other parameters, (c) use of 
direct bearing and elevation angle observations in combination 


with observations of radio frequency (RF) energy reflected 


LO 





from the water or Earth's surface, and/or (da) a series of 
observer maneuvers used in conjunction with angular observations. 

Determination of a target's position and motion by means 
of calculations performed on observations of only bearing 
(and possibly elevation) angles is confounded by the reality 
that the angle measurements are noisy. That is, while a 
series of measurements may be centered on the true values, 
individual measurements will differ from the true values by 
varying amounts and in a random way. This type of PR prob- 
lem may be thought of, then, as having two major facets: 
first, from the series of measurements the observer must 
make "best guesses" as to the true angle values; and second, 
using these estimates, he must solve a dynamic geometry 
problem to ascertain target position and motion. 

A "good" system for passive ranging will determine 
target parameters using only its own observations, but it 
might reasonably be expected to also make use of any a 
priori targeting information to improve its speed and/or 
accuracy of solution. A method for doing this is by means 
Sena Kalman Filter (KF) algorithm in which the prior infor- 
mation is utilized along with an estimate of how accurate 
that information is. The KF (or, more precisely, the Ex- 
tended Kalman Filter (EKF)) is successfully employed presently 
in the submarine and surface ship anti-submarine warfare 
communities. In this environment observations are normally 


acoustic bearing measurements, although the use of 
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depression/felevation angle with "bottom-bounce" reflections 
has been utilized or proposed. 

Representative of current thinking on passive ranging as 
an air-to-air problem are techniques proposed by Ohlmeyer 
[Ref. 1], and Stiffler [Ref. 2], which also rely on the use 
of EKF's. Ohlmeyer's technique uses only bearing angle 
information, whereas Stiffler presents six different algorithms 
using bearing-only or a combination of bearing and elevation 
angle observations. 

This paper proposes and investigates methods of passive 
target motion analysis by means of two EKF's, where both 
bearing-only and bearing-elevation solutions are computed 


Simultaneously for comparison purposes. 


I 





II. THE KALMAN FILTER AND EXTENDED KALMAN FILTER 


ieee tHE KALMAN FILTER 

The following discussion briefly describes Kalman Filter- 
ing in general, the need in our case for the Extended Kalman 
Filter, and an application of the EKF to the air-to-air 
passive ranging problem. The reader interested in a more 
in-depth development of the KF and EKF is referred to 
References 3, 4, and 5. 

The purpose of a KF is to keep track of the state of a 
system by means of a sequence of noisy measurements. This 
is done by recursively updating an estimate of the state by 
computationally processing the succession of noisy observa- 
tions in such a way as to reduce as much as possible the 
errect of measurement errors. For our purposes we will 
assume that the object (target) whose state we are esti- 
mating, and the observer (fighter) exist in a Cartesian 
coordinate system and that the state, S, is composed of a 


four- or five-element vector of the form: 


or S = (21) 


KM Ke xe 
PM NM Ke Xe 


fe 





where X and Y are the X and Y components of target velocity, 
X, Y and A are the X, Y and altitude components of target 
position; and where the true values of X, Y and A are assumed 
constant. Note that the KF in general is in no way limited 
to this type of reference frame or state vector, and that 
observations from which the estimated state is computed need 
not have the same dimensionality or units as the state 
vector. In our case, for instance, each observation will be 
either a bearing angle measurement or a vector of both 
bearing and elevation angle measurements. 

The KF is a predictor-corrector type estimator that 
propagates an estimate, u, of the state, S, and an associated 
covariance matrix, ), reflective of the accuracy of the 
state estimate. The KF algorithm is implemented in two 
alternating stages. First, previous estimates of u and } 
are extrapolated one time-step ahead based on the assumed 
system dynamics; this we will refer to as the Movement Step. 
The extrapolated quantities are then used to compute a set 
of optimum weights called Kalman Gains (simply proportion- 
ality factors). These are applied to the prediction and to 
a new observation in a Measurement Step, which provides an 
updated (corrected) estimate of the state and its covariance. 
The process is then repeated. 

Symbolically, the KF algorithm may be written as follows: 


Mevement Step: 


uw. = Ou, tree tl ZZ) 





= ¢jy. , o + Q (ES) 


J j-1 
Kalman Gain computation: 
K, = ; a” UE , fee P) (2.4) 
Measurement Step: 
Ms = ou, + Ki(Z~ uy ~ Hus) (2.5) 
7 = KH]. (Za0,) 


The notation is as follows: 

@’ is the transition matrix which describes how the 
System state changes between observations, 

H is the transition matrix which describes how the 
measurement depends on the state, 


uo, Q, u_, and P are respectively, movement noise mean 


Ww V 


and covariance, and measurement noise mean and covariance. 
(Note that movement noise and measurement noise are assumed 
to be normally distributed random variables, i.e., 
W ~ N(u..Q) and V ~ N(u,P)). We will generally assume here 
that Met; Q, and WL are zero. 

K is the Kalman Gain, 

Zis the observation (measurement), 


I is the identity matrix, 
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J is the observation time index, 
T (aS superscript) indicates transpose, 
~ and « indicate estimates after Movement and Measure- 


ment, respectively. 


B. THE EXTENDED KALMAN FILTER 

The above described KF is appropriate when both the 
movement and measurement models are linear functions of the 
state variables. We have, however, a non-linear measurement 


Bunction since 


Z = h(S) + V 


where the measurement, Z, is a function, h(S), of system 
state plus some noise (error), V, and where h is composed of 
one (or two) arctangent functions of state vector elements. 


We have: 


h(s) = tan”? + | in the bearing-only case, and 
AX 
tan AY 
h(S) = in the bearing- 
tan au 
(ax*+ay*) 272 J elevation case, 


where AX, AY, and AA are, respectively, the X, Y and altitude 
differences between fighter and target positions. 
Since the measurement in both the bearing-only and 


bearing-elevation cases iS a non-linear function of the 
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Current state vector, the measurement transition matrix H 
in the KF algorithm must somehow be "linearized". This is 
accomplished by expanding h in a Taylor Series about a known 
vector that is close to the true state S. For this known 
vector we will use wu, our current estimate of S. Truncating 


after the first two terms of the Taylor Series we have that 


tS ee tl + cles) 7 oS 
SS 
The effect on the KF algorithm is that in Equation 2.5 
ie is replaced by the function Athy.) which is the measure- 
ment vector as a function of the state estimate, and in 
Mmem@ation 2.4 and 2.6 H = 0h(S)/dS . One condition re- 


S=yu 
quired here is that the estimate be unbiased, i.e., u_, 


Vv 
the estimated value of V, must be zero [Ref. 3]. 

This extension of the KF to the non-linear case is 
appropriately called the Extended Kalman Filter (EKF). 
Higher order, more precise approximations to the optimal 
non-linear filter can be achieved using more terms of the 
Taylor Series expansion for the nonlinearities, and by 
deriving recursive relations for the higher moments of S. 
See Reference 3 for a discussion of methods of this type. 

The (non-extended) Kalman Filter is optimal when its 
assumptions are valid [Ref. 6], however, this is not neces- 


Sarily true of the EKF. Since, in our case for instance, 


the measurement transition matrix H both depends on the 


ae, 





current state estimate and is used to obtain revised state 
estimates, there is a potential for bad estimates to get 
worse and for complete loss of track to occur [{Ref. 7]. 
Also, the matrix ) is only an approximation to the true 
covariance matrix and is not always an accurate indicator of 
the estimation errors. Fortunately, however, the EKF has 
been found to perform well in a number of practical applica- 
tions [{Ref. 3] including passive ranging. 

The design of practical, operational filters "is an art 
wherein one attempts to retain the behaviour of the EKF when 
things are going well, while simultaneously being able to 
recognize and correct for incipient loss of track" [Ref. 7]. 
Techniques for doing this (and cautions) are discussed later 


and in References l, 3 and 8. 


C. FUNCTIONS AND MATRICES 

Functions and matrices are given here for the bearing- 
elevation cases. Those for the bearing-only cases are 
exactly the same, but excluding the last (fifth) element 
in the case of a vector, and the last row and column in the 
case of a matrix. 

With u as the estimate of the state vector 


o— ow ,.X, val, we have that the movement function 
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a i-l 
Yi ~ ea 
g(S) Rg <P Sy 5) ammllie 
S=u 
eee eee AG 
2H eee 5 
so the movement matrix is 
1 0 0 0 0 
0 ik 0 0 0 
6 = iL 0 iL 0 0 
0 i 0 iL 0 
0 0 0 0 iL 
The covariance matrix is 
g2 0 0 0 0 
a: 2 
0) Oy 0 0) 0 
2 
) _ 0 0 Oy 0 0 
2 
0 0) 0 oy 0 
2 
0 0 0 0 on 
The measurement function is 
one ieee) 
tan 1 cae 
ia ora: 
hits) = 
Sa ee ea 


Up | 


iby) 





a a _ als 2 
R > (Cu, F..) (Hy FY) ) , and where Poy Fs 


F. are, respectively, the fighter's known X, Y, and altitude 


where yu and 


positions, so the measurement matrix is 





-F UL oe 
0 oO ve = * = 0 
jal u 
H(S) = e . 
S=u ; . = (Ux-F,) (Hy-FY) —(u,-F) (UF) Up 
a ee ee ne) 
(USt(Ua Fa) Yup  (ueH(uy-Fy) Duy uot (uaF) 
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IIIT. PARAMETERIZATION AND ASSUMPTIONS 


A. FIGHTER AND TARGET FLIGHT PATHS 

For all simulated tracking sessions the fighter utilizes 
a "Zig-zag" pattern as described in Appendix B and illus- 
trated in Figure 3.1. Briefly, this pattern has the fighter 
initially turn so as to have the first bearing observation 
0.7 radians (about 40°) right of dead ahead, hold this 
heading for 1 minute, then turn to put the most recent bearing 
0.7 radians left, now continuing for 2 minutes. Alternating 
2-minute legs are continued until a specified number of 
observations have been taken, or until the target passes by 
the fighter. The target, meanwhile, in all cases maintains 
a non-turning constant altitude track. For the computer 
Simulation runs each tracking session is limited to 5 minutes 


(61 observations at 5-second intervals). 


Bs. INITIAL PARAMETERS 

Eighteen different combinations of P and ) values (that 
is, covariances for observation error and state estimate) 
are utilized for initial parameterization in the simulations. 
The runs are first done using the "flat-Earth" motion model, 
then they are all repeated using the spherical-Earth model 
with an elevation angle correction (Appendix D), then twelve 


cases are simulated using the spherical~-Earth model without 
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the elevation correction, and finally, three cases are simu- 
lated where a 5° constant bias is added to the angle measure- 
ments. 1000 replications of each set are run for both the 
bearing-only and bearing-elevation cases. In all cases the 
target, at time 0, is at 80 nautical miles (nm) range, heading 
directly toward the fighter at an altitude of 7nm (42,532 
feet) and a ground speed of 500 knots. The fighter maintains 
an altitude of 3nm and a ground speed of 300 knots. The EKF 
initial estimates of target state for each run are varied 
randomly according to a multivariate normal distribution 
centered at the correct values, and with covariance as stated 
below. 

At each observation the bearing-only (BO) and bearing- 
elevation (BE) EKF's are fed identical random azimuth 
measurements, while the BE EKF additionally gets random ele- 
vation angle measurements. The random measurements are 
generated by adding to the true angle a random number from 
a normal distribution with mean 0 and variance P. 

The P-values used are —- Qe. La ae ie, and foo 
degrees“, where these values apply to both bearing and, 


where applicable, elevation measurements. For example, in 


em@e first set of bearing-only simulations P = [25] whereas for 
the bearing-elevation runs P = WS elke We assume here 


that both bearing and elevation measurements have the same 
Megiteude and distribution of errors (normal), and that 


they are independent of each other. 
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The initial estimate of the covariance matrix has diagonal 


values: of = 1000 hao BSS oe = 250,000 knots“, of = Anm-, 
and x = Ara , while a is varied from 1600nm* to MOO aC 


A a 


1e. 100nm7 for the different cases. The choices for these 
values as well as for the P-values are somewhat arbitrary 

and their validity will not be argued here. The reader is 
referred to Appendix B for a description of how the reference 
frames (coordinate systems) are defined, and of how true 


positions and observation angles are computed. 


C. ASSUMPTIONS 
For EKF computer simulations the following assumptions 
are made: 

(1) The observer (fighter) knows his own position, 
course and speed precisely. 

(2) Maneuvers are instantaneous; there is neither time 
loss nor positional change associated with the 
observer's course changes. 

(3) The target holds a constant course, speed and 
altitude, while the fighter holds constant speed 
and altitude. 

(4) There is no refraction or other distortion of the 
radiated energy. 

(5) Angle measurement errors in elevation and in azimuth 
have a Gaussian (normal) distribution with zero mean 
and known variance. Errors are assumed to be inde- 


pendent of range. 
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IV. SIMULATION RESULTS AND OBSERVATIONS 


A. TABLES 

Results of the computer simulations are summarized in 
Tables I-IV which are for, respectively, the planar ("flat- 
Earth") motion model, the spherical-Earth motion model with 
elevation angle correction, the spherical model without 
elevation angle correction, and, finally, the spherical 
model with elevation angle correction, but with a constant 
angle bias inserted into the observations. Data in the 
tables are categorized by initial range variance, and for 
Gach range variance the different selections of angle measure- 
ment precision are represented, with data from the bearing- 
@maly EKF on the left and from the bearing-~elevation EKF on 
the right. 

Comparison data are the proportions of dropped tracks 
using three different criteria. First is the "e€ > 30" column, 
meen lists the proportion of runs in which actual range 
error was greater than the 3 standard deviations as estimated 
by the ERKF, 1.e., € = |R-R| > 30. (See Appendix A for a 


description of computations of R and om ) The "Rg" column 


R° 
lists as dropped tracks those for which the final range esti- 
Mate 1s in error by more than 25% of the true range, and the 


"Rg-Vvel" column lists those for which either range or 


velocity are in error by more than 25% of actual values. 
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TABLE i 
Simulation Results--"Flat-Earth" Motion Model 
Proportion of Dropped Tracks 
(by criterion) 


1000 Replications 


Bearing-only Bearing & Elevation 
Initial Range Variance = }(4,4) = 40°nm* 


1g Degrees 


Measurement 
- O Re cul Rg Precision E>30 Rg-Vel Rg 
0.188 0.949 0.732 ' 0.075 0.880 0.554 
moss 0.757 0.452 2 0.047 0.473 0.169 
0.094 0.480 0.175 ib 0.073 0.206 0.059 
Ma24 0.235 0.109 Se Gime orisg) 0.076 
M17 4620. 179~Ss«éO. 128 ie 0.208 0.153 0.092 
mmee7 )60OCO0.177~S—«0. 138 02 Omeao er Geies | 6.006 
y (4,4) = 20°nm? 
ZG} Spa IL Rg P E>30 Rg-Vel Rg 
0.121 0.904 0.710 5 0.042 0.781 0.504 
Mmol7 0.664 0.329 2 0.004 0.347 0.096 
9.006 0.328 0.050 iy 0.004 0.038 0.001 
mola 0.083 0.001 35 0.006 0.009 0.001 
0.066 0.029 0.007 es 0.041 0.016 0.003 
0.241 0.029 0.007 02 0.184 0.012 0.001 
¥ (4,4) = 10°nm? 
a oO Rg-Vel Rg 2 Se sio, ine Wretil Rg 
Mes, 0.858 0.713 55 0.032 0.701 0.484 
Mol? 0.502 0.293 2 0.002 0.203 0.071 
0.011 0.239 0.030 Ik 0.001 0.018 0.000 
0.003 0.039 0.000 55 0.004 0.000 0.000 
0.014 0.002 0.000 ie 0.006 0.000 0.000 
melo? 0.001 0.000 02 0.037 0.000 0.000 
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B. OBSERVATIONS 
Some immediate observations can be made from the tables: 

1. From Table I ("flat-Earth" Motion Model), the bearing- 
elevation results appear superior to those for bearing-only 
for every set of parameters when using either the range- 
velocity or range error criteria, and they are generally 
superior when using the € > 30 criterion as well. It should 
be noted, however, that in all of the cases represented the 
fighter-target altitude difference is held at 4nm (about 24,000 
ft). The advantage of BE over BO was shown by limited addi- 
tional simulations to disappear as the altitude difference 
mame CO Zero. 

2. Improvement in the accuracy of the initial range 


estimate (decrease in the variance ray) fcem 407 ae 20, 


and then to 10°nm- has a positive effect in decreasing the 

proportion of dropped tracks as one might reasonably expect. 
3. Improvement in angle measurement precision (variance 

P) generally improves the EKF's ability to track. It seems 

evident, for instance, that P = 5° is impractical from an 

Operational standpoint since neither EKF could determine 

target range and velocity within 25% even 30% of the time 

regardless of initial range variance. It is interesting to 

note that the ¢ > 30 dropped track criterion for the BE 

mee for P = 3c is exceeded in a much lower (2-10%) proportion 


of the cases than are the range and range-velocity standards 


Mee-J5%). This is a valuable trait in that, although the 
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pearing-only 


E>30 


E255 
0.082 
0.076 
O25 
OE 231 
0) 


e>30 


0.238 
0.024 
0.004 
0.003 
0.012 
0.088 


TABLE If 


Simulation Results--Spherical Motion Model 
with Elevation Correction 


PeepornelOmsOnee:opped Tracks 


(by criterion) 


1000 Replications 


Initial Range Variance = } (4,4) = 40*nm* 
P Degrees” 
Measurement 
Rg-Vel Rg Precision E>30 Rg-Vel 
0.937 0.751 55 0.121 0.878 
0.718 0.370 2 0.073 0.466 
0.458 0.145 ibs 0.085 0.159 
0.224 0.106 & 0.160 0.139 
0.178 0.120 Elis 0.440 0.179 
0.170 0.129 02 0.858 Oneee 
Y (4,4) = 20°nm? 
Rg-Vel Rg P E>30 Rg-Vel 
0.919 0.752 x 0.122 0.808 
D,G28 0) aap 2 0.010 0.291 
0.350 0.049 iL, 0.012 0.043 
0.084 0.007 5s 0.034 0.013 
0.027 0.008 sips 0.247 0.012 
0.028 0.006 Oe 0.801 0.078 
(4,4) = 10°nm? 
Rg-Vel Rg P E>30 Rg-Vel 
0.842 0.721 56 0.083 0.663 
Gems) 01232 2 0.003 0.203 
DPI ove ibs 0.005 0.015 
0.047 0.000 = 0.006 0.000 
0.001 0.000 ce 0.142 0.000 
0.000 0.000 02 0.742 0.059 
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Bearing & Elevation 


Rg 


ORS9 2 
5195 
0.049 
02075 
Oyen 
0.174 


0.542 
ORO 
0.004 
07003 
0.006 
0.045 


0.478 
0.066 
0.000 
0.000 
0.000 
0.050 


EKF estimates of target state are not very accurate, the 

EKF iS apparently "aware" of this since its estimated range 

variance is evidently high. Since the state covariance 

matrix iS available to the user, he should generally be 

able to determine the worth of the estimated state information. 
4. This characteristic does not necessarily hold true 

when measurement precision is extremely high, however; in 

fact, it becomes less true as P-values are decreased below 
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for these EKF's (e.g., Table I, )(4,4) = 40° and P = .02%). 
A similar observation was made in Reference 6. This apparent 
EKF trait may be explained by the Taylor Series approxima- 
tion used to linearize the Measurement Function (Chapter 

II) inducing more error than that added by the noisiness of 
the bearing measurements. The EKF-estimated range variance 
becomes relatively small because of the precise measurements, 
Bue the actual range errors are not correspondingly reduced 
apparently because of the EKF linearization approximation. 
Practical results are still quite good as evidenced by the 
still relatively small proportion of dropped tracks using 

meen range and range-velocity criteria. With a longer obser- 
vation period, however, we might expect the estimated 

range variance to continue its decrease relative to the 
Squared range error (covariance collapse), and we would 


anticipate more tracks being dropped. A possible solution to 


this problem, if highly precise angle measurements are 
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TABLES LTT 


Simulation Results--Spherical Motion Model 
Without Elevation Correction 


Proportion of Dropped Tracks 


(by criterion) 


1000 Replications 


Bearing-only 


Bearing & Elevation 


2 


Initial Range Variance = }) (4,4) 40“nm 
P Degrees 
Measurement 
E>30 Rg-Vel Rg Precision E>30 Rg-Vel Rg 
0.079 0.726 8G 2° 0.079 0.477 0.126 
0.071 0.420 0.127 1 0.137 0.239 0.055 
0.108 0.205 0.090 5° 0.372 0.130 0.071 
0.219 0.017 0.129 ae 0.998 0.166 0.117 
¥ (4,4) = 20°nm? 
E>30 Rg-Vel Rg P E>30 Rg-Vel Rg 
0.035 0.609 0.254 2° 0.028 0.336 0.044 
0.069 0.156 0.002 1 0.069 0.156 0.002 
0.011 0.073 0.006 ee 0.312 0.022 0.002 
0.053 0.017 0.002 Ae 1.000 0.005 0.002 
yea) = 107nm* 
E>30 Rg-Vel Rg P E>30 Rg-Vel Rg 
0.023 0.502 0.237 2° 0.027 0.226 0.050 
0.002 0.222 0.017 1 0.065 0.098 0.000 
0.006 0.041 0.000 ee 0.328 0.019 0.000 
0.008 0.001 0.000 ac 1.000 0.000 0.000 
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TABLE IV 


Simulation Results--Spherical Motion Model 
with Measurement Bias 


EqO UOrLICneOLeDroppedatracks 


Jesu iaie iia Sign copgy 


1000 Replications 


Bearing-only Bearing & Elevation 
Initial Range VAriance = ) (4,4) = 40°nm* 
Measurement Precision = P = 54 

AZ Elev 
E>30 Rg-Vel Rg Bias Bias E>30 Rg-Vel Rg 
ess 0.238 0.113 0 5 Oz220 0.838 ORZS5 
Mees! 60.999 0.989 3) 0 0.865 O29ou 02752 
Meeo> 0.999 0.986 5 5 0.747 02298 0.749 


available, would be to use a higher order Taylor ties 
approximation as previously suggested. 

5. A comparison of Tables I and II indicates that spheri- 
cal motion target data can be accommodated by the flat- 
Earth EKF when an appropriate elevation angle correction 
1S applied. At the very high precision levels, however, 
there 1s a notable degradation with the BE EKF results, 
especially when using the ¢€ > 30 criterion. This is apparently 
@ue tO the filter being “confused” by highly accurate 
measurements which are "corrected" by an elevation angle 


correction (for Earth's curvature) which is no longer 


3. 





adequate. The result is that estimated range variance is 
considerably smaller than actual range variance. 

6. A comparison of Tables II and III indicates that the 
BE EKF is somewhat improved when using the elevation angle 
MBeerection to allow for the Earth's curvature, but the filter 
without the correction still performs almost as effectively. 
(The BO EKF, of course, is not affected by the correction 
Since it does not use elevation angle as an input). As 
angle measurement precision increases, the positive effect 
of an elevation angle correction to the BE EKF also increases, 


however (e.g., Tables I and II, ) (4,4) = 40%, P=: oi 


Ne 

7. Table IV illustrates the effects of biased measure- 
ment data on EKF performance. When azimuth angle, or both 
azimuth and elevation angles have a bias, then neither 
filter operates effectively. When only the elevation angle 
1s biased the BE EKF performance 1s severely degraded. Only 
the P = .25 measurement precision was investigated; it was 


felt to be a good test case since it consistently produced 


among the best results with unbiased inputs. 


eee Ll LME SEQUENCE PLOTS 

Figure 4.1 is a set of plots showing a time sequence of 
range errors during the 5 minutes of tracking. The median 
values of range errors are represented for the 1000 repli- 
cations of each of the different levels of measurement pre- 


cision (P). A representative set of plots are shown in 
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Figure 4.1. Median Range Error Vs. Time--"Flat-Earth" 
Model, Range Variance = 100nm 


So 


miciGomowmewitch are Ger the "“flat—-Earth" motion model with 
range variance of 100nm*. Similar plots for the other sets 
of parameters are included at the end of the chapter. Note 
that the horizontal axis is not always at zero range error, 
and that the vertical scales differ among plots. In each 
plot the solid line represents BO EKF results, whereas the 
dotted line represents those from the BE EKF. 

Several important features of the EKF's are illustrated 
by the plots: 

l. The BE EKF is again seen to be generally superior to 
the BO EKF. The major advantage of the BE EKF seems to be 
the generally greater speed with which the estimated range 
converges to true range particularly in the lower measure- 
ment precision cases. At very fine measurement precision 
(P = .01 and P = .0004 degrees“) there seems to be little 
if any advantage to using both bearing and elevation angles. 

2. In all cases the trend seems to be toward better 
convergence with time, once the initial transients have 
disappeared. This suggests that even in those cases where 
a large proportion of the runs are "dropped tracks" according 
to the range or range-velocity criterion, additional tracking 


time, if available would bring increasing numbers into 


a 
convergence. Limited additional computer simulations seemed 
to bear this out. 

3. The “humps" which are especially prominent in the 


a ae i a nil 5° cases occur at about 2 and 4 minutes. 
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They indicate that range error is decreasing from time 0 to 
1, increasing from time 1 to 2, decreasing again from 2 
Mons, etc. Referring to Figure 3.1 (Fighter and Target 
Flight Paths) the convergence of estimated range to true 
range, then, appears to increase as the rate of bearing 
change increases. This is consistent with the well known 
(among submariners) result that in passive ranging "the 
greater the time rate of change of bearing the more accurate 
the range estimate" [Ref. 9]. This rate is slowest (with 
this particular geometry) when the fighter first turns in 
toward target flight path, increasing until he next turns in 
from the opposite side of target flight path. The humps, 
then, roughly coincide with the times at which the fighter 
1s turning back in, 1.e., the times when bearing change 

rate is relatively low. 

4. It is evident that the time to convergence is greatly 
affected by angle measurement precision. Note that in the 
P = Ro case convergence is extremely rapid after 1 minute, 
which is the time of the first fighter turn. This suggests 
that in the case of extremely good precision, convergence 
time can be even further decreased by making the fighter's 
first turn earlier, or perhaps having the fighter maintain 
a constantly turning flight path such as the "Sine wave" 
flight path proposed by Ohlmeyer [Ref. 1}. Limited addi- 
tional simulations with the fine precision and an earlier 
turn did, in fact, indicate that convergence time would be 


reduced. 
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Figure 4.5. Median Range Error Vs. Time--Spherical 
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Figure 4.8. Median Range Error Vs. Time--Spherical 
Model (No Elev. Correction), Range 
Variance = 400nm 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Computer simulations using the two EKF algorithms 
presented here demonstrate that both have the potential 
to be valuable tools for air-to-air passive ranging as 
long as their assumptions are met. The bearing-elevation 
EKF generally performs better than the bearing-only EKF 
when the fighter and target have a significant separation 


in altitude. 


B. RECOMMENDATIONS 

There are several possible EKF modifications which seem 
to offer potential in improving its performance. Some of 
these are: 

l. Filter re-initialization when it is determined to 
be grossly in error. Aidala [Ref. 8], discusses ways in 
which "premature covariance collapse" can be recognized. 
He notes that onset of this malady, which usually leads 
to a dropped track, often occurs after only two initial 
measurements. This suggests that if it can be recognized 
there would normally be ample time to restart the EKF. He 
also cautions that some "previous attempts to prevent ab- 
normal filter behavior have proved ineffectual when sub- 


jected to rigorous testing...because they focus only on 
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the symptoms of the problem (rather than) the underlying 
causes" [Ref. 8]. 

2. Filter self-correction when obtaining "unbelievable" 
results. Although this may fall into the category of 
treating the "symptom" rather than the ailment, it seems 
reasonable to put some bounds on the acceptable estimated 
target state based on physical knowledge. The most obvious 
limit would apply to a low target; for example a fighter 
at an altitude of 3nm with a target elevation observation of 
-4° could immediately deduce that the target is (depending 
on measurement precision) no farther away from him than 
Bees CSC 4° = 43nm, since he can generally be assumed to 
be above the Earth's surface, so it would seem futile to 
waste more time estimating him at, say, 80nm. Similarly, 
the maximum altitude of a target may be known, as well as 
possible limits on speeds, ranges, or headings. These 
things would be scenario-dependent and would probably re- 
quire aircrew inputs and reasonable default values. 

3. Increased frequency of observations (filter cycle 
rate). A five-second interval was chosen for these simula- 
tions because that seemed long enough to ensure independence 
among measurements in an operational system, and because it 
reduced the computer time required to do the many replica- 
tions. A few simulations using three-second and one-second 
intervals obtained slightly faster convergence though little 


improvement in the numbers of dropped tracks. 
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4. Optimized fighter PR maneuver. Since the rate of con- 
vergence to the true target state evidently increases with 
the rate of bearing angle change, this should be considered 
in choosing fighter maneuvers. A pure collision course is 
not conducive to good passive ranging, especially for the 
BO EKF, since it precipitates no bearing change, but an 
alternating lead-lag, for instance, may accomplish virtually 
the same thing tactically, yet allow sufficient bearing 
change for a PR solution. The "optimum" maneuver would 
undoubtedly be different for varying scenarios and equip- 
ments. Altitude and speed changes as well as horizontal 
turns should be considered. The APL program (Appendix E) 
used to model these EKF's can be easily modified to accommo- 
date different fighter and target flight paths, and so could 
be a useful tool in comparing PR maneuvers. 

5. P-values (observation precision) adjusted to actual 
equipment specifications. This may entail varying the 
P-values as a function of range. Simulation results using 
this change are shown in Reference l. Other equipment-peculiar 
Characteristics such as own-ship navigation error would also 
have to be considered. For these simulations it was assumed 
that the distribution of angle measurement errors was known. 
The effects of giving the EKF wrong precision error distri- 
bution information should be investigated. 

6. Interactive system. The aircrewman may have a 


variety of sources of a priori information with which to 
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make an initial or updated estimate of target state. There 
should be a means for him to insert this information along 
with a confidence level or estimated variance, and, as 
mentioned above, there should be a set of reasonable default 
values. 

7. Measurement variance reduction. For example, several 
measurements could perhaps be averaged to obtain one obser- 
vation input which is less affected by individual outliers. 
This could be particularly important for the first few 
observations, since as Aidala [Ref. 8] suggests, they are 
so influential in the "premature covariance collapse" 
problem. 

8. Maneuvering target recognition and tracking. This 
continues to be an unsolved problem. Ohlmeyer [Ref. 1], 
Suggests "that small levels of target acceleration can be 
accommodated within the filter (uSing methods such as) the 
introduction of state process noise.” 

9. Parallel filters. If computer capacity is available, 
it may be advantageous to have more than one EKF operating 
at the same time on the same target. The separate filters 
Mmegnt use alternating observation inputs, or perhaps initialize 
at different times. Besides rendering a "Second opinion" 
this could be helpful in the maneuvering target situation. 

10. Multiple target processing. Since fighter maneuvers 
need not generally be precisely designed for a specific 


target, the fighter could be making his maneuvers against 
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@ particular, perhaps primary, target, and concurrently be 
computing passive ranging solutions on other targets if 

his equipment has sufficient sensor capability and computer 
Capacity. 

li. Elevation angle-only EKF. There could be situations 
mer which the fighter would want to maintain a pomecron 
course in azimuth (precluding bearings-only passive ranging) 
and yet desire a PR solution. This could conceivably be 
accomplished with altitude changes and an EKF which requires 


only elevation angle observations. 
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APPENDIX A 


VALIDATION 


A. VALIDATION METHOD 

The method of validation used here is borrowed from 
Daly [Ref. 6]. In essence, validity is claimed if the EKF 
tracks in range variance as well as in range. 

If the EKF is operating properly andifit is being 
provided with accurate enough information, then the esti- 
mates of range variance after individual tracking sessions 
should be very close to the true range variance. Compari- 
sons of the EKF variance estimates over many tracking 
sessions with the actual (sample) variances provide a means 
of evaluating the consistency of the EKF range estimates. 
As a measure of this consistency we will use the statistic 
Kl which is the ratio of sample range variance to EKF- 
estimated range variance at the end of a tracking session, 


averaged over n distinct sessions. 


Pas 


R, is the EKF range estimate, R. is the true range, and 


of is the EKF range variance estimate, all at the end of 
0 
tracking session i. The tracking session is simulated 
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n = 1000 times using the same values for all input param- 
eters, but using distinct random number streams for measure- 
ment errors and initial state estimate errors. The method 
Gor Computing R and oe is described below in Section C. 

We expect that if our theory is correct and the EKF 
has accurately estimated the true range variance, that the 
value of Kl will tend toward unity. 

The EKF algorithms were applied to the same passive 
ranging problem (scenario) which was addressed in Reference 
6. In this problem the observer and target are both surface 


ships. At time 0 the target bears 090° at 20nm. The target 


proceeds on a course of 360° at 10 kts, while the observer, 


TABLE V 


Kl Values--Surface Ship Scenario 


Flat-Earth Motion Model 


1000 Replications 


2 


P idea) =? 10 ea) = 7 
Ns 0.982 1.015 
ee 0.957 0.015 
Lo 0.991 0.953 
ome 0.985 0.959 
ome 0.947 0.976 
On 1.048 0.955 
none Lalas 1.091 
mon 1.140 1.086 
Lone 3.750 0.941 
omme pon Ee: 
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for his passive ranging maneuver, proceeds on a course of 
180° for 15 minutes, thenturns (instantaneously) to a course 
of 360° for another 15 minutes. A total of 11 observations 


are taken at 3 minute intervals. 


Be VALIDATION RESULTS 

Table V illustrates that the Kl values do indeed 
generally tend toward unity, thus supporting our claim of 
validity. It is notable that with extremely fine measure- 
ment precision the Kl values no longer stay near 1. This 
finding was also reported in Reference 6. As discussed in 
Chapter IV, this seems to be a result of the range variance 
decreasing because of the fine precision, whereas the actual 
range errors are apparently limited in their decrease by 
the errors inherent in the EKF due to its linearization 
approximation in the Measurement Step. 

Results in Reference 6 were obtained using an EKF 
parameterized with range, bearing, velocity, and direction, 
faemer than X, Y, X-velocity, and Y-velocity as used here. 
The author reports that he had obtained similar results 
using the latter parameterization, but that when using the 
(X,Y,X,Y) parameterization the Kl values stayed close to l 
for a higher degree of measurement precision. Results con- 
tained here appear to confirm that. 

When this same technique was applied to the scenario 
described in Chapter III the resulting Kl values were less 


clearly convergent to 1 (Table VI). While being somewhat 
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TABLE VI 


Kl Values--Fighter Scenario 


Flat-Earth Motion Model 


1000 Replications 


Bearing-Only 


he Mm FR NY WM 


cantly more than those for the simpler problem 


P 5 (4,4) = 40° ¥(4,4) = 20° 5 (4,4) = 10° 
5 0.981 0.893 0.897 
2 0.932 0.882 0.873 
1 0.904 0.868 0.835 
5 0.853 0.840 0.833 
1 0.851 0.799 0.834 
02° 0.938 0.843 0.826 
Bearing & Elevation 
P 5 (4,4) = 40° ) (4,4) = 207 5 (4,4) = 10° 
2 1.090 0.964 0.900 
2 0.893 0.846 0.832 
é 0.839 0.826 0.837 
2 0.831 0.801 0.774 
2 0.889 0.816 0.823 
02° 1.146 0.983 0.829 
"close" to 1 (+0.2) the Kl values varied from 1 signifi- 


Gaal Pe Os) ee 


Only those values which were less than 10 were included in 


this computation of the average Kl's; this was done to 


eliminate the effect of those 


(normally few) dropped track 


6 


Kl values which diverged drastically from 1, e.g., 10 , 
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due to the much larger number of iterations per tracking 
session. This limit was picked arbitrarily and may partially 
account for most of the average Kl's being less than l. 

The variation in the values seems to infer that the rate of 
convergence is scenario-dependent in that it is affected 


by the speeds, distances, and geometry involved. 


C. COMPUTATION OF RANGE VARIANCE 
The covariance matrix for this type of parameterization 
has terms for X, Y, and A (as well as X and Y) but not for 


range itself, so the variance of the range, of 1s not 


R?” 

provided directly, but is obtained as follows: 
For the bearing-only cases, with R, Ry, and Ry repre- 

senting range, and X-and Y-components of range respectively 


_ 2 ME//2: 
R = (Ry + RY) 


Q 


R 
SO = =, and aR x, 


dRy 
The first-order (linear) Taylor Series approximation 


(where ~*~ indicates an estimated value) is: 


cas 


re oo oR es OR ne 
R = R + oye ae R,,) + OR, Fy Ry) 
therefore 


2 E oR .2 2 oR oR oR 2. 2ZeZz 
e} = (———) “go. + (——) (~—) Cov ( ; ) + (—) oa 
R OR Xx aR aR, Ry Y aR, Y 


ys 





or, substituting for the partial derivatives, 


a eee 2 2 
on = BEAK x + 2R,R,Cov(R,,R,) + Roy): 


Similarly, for the bearing-elevation cases, since 


R = (RE + RS - seg? 


then 
2 ae epee = 2 2a? 
on = BRC ar yy ye yea Rye Ry) 
~ R-o“+2R ReeGor (ines) ) 
enya x a Ge eae 
Note that Cov (Ry, Ry) = Cov(X,Y) where X and Y are the target's 


X and Y coordinate positions, since the fighter (observer) 


position is known. 
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APPENDIX B 


EKF REFERENCE FRAMES AND MOTION MODELS 


A. "FPLAT-EARTH" REFERENCE FRAME 

The EKF’s for both the bearing-only and bearing- 
elevation observation models are set in a Cartesian 
coordinate system with axes X, Y and A. This reference 
system has its origin at 0 altitude directly beneath the 
fighter at time 0 (time of initial observation), and Y-axis 
in the direction of the first observed bearing. Bearing 
and course are measured in the X-Y plane clockwise from the 


Y-axis. 


B. MOVEMENT IN THE "FLAT-EARTH" MODEL 

Since altitude jis held constant for the fighter and the 
target, movement is a two-dimensional problem. For example 
if a target is at position (X,Y) A) at time 1, and at time 
2 has moved D units on a course of C, then his new position 


(X5,Y.,A,) is computed by: X, = X, + D sinc, Y, = Y 


1 1 + Dcosc, 


and A, = A 


5° 
C. SPHERICAL REFERENCE FRAME 

Although the EKF algorithms presented here utilize a 
Cartesian (flat-Earth) reference frame, the algorithms are 


also tested using fighter andtarget positions generated by 


a spherical model. Here the Earth is modeled as a perfect 
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sphere of radius 3444 nautical miles. "Constant" course in 
this model means a great circle track with constant alti- 
tude; speed refers to ground speed, i.e., the speed of the 
object's projection on the Earth's surface. 

For computational ease a reference system is again based 
on fighter location and bearing to the target at time 0. 
The origin of this (latitude, longitude, altitude) spherical 
coordinate system is defined as the zero-altitude point 
directly beneath the fighter at time 0, and the initial 
observed bearing to the target is defined to be 0. (This 
may be thought of as shifting the Earth's coordinate system 
so that the intersection of the equator and the Greenwich 
Meridian is directly beneath the fighter at time 0, and the 


first observed bearing is True North). 


D. MOVEMENT IN THE SPHERICAL MODEL 
In computing movement in the spherical system the 
following rules of spherical trigonometry are used (refer 


mompagure B.1): 


Law of Sines: 


Solem. _ Sin 8 Lee San. ¥ 


Sin AJ/E  sinB/E sin C/E 
Law of Cosines: 
2650578) =) COs A/JE Cos C/E + sin A/E sin C/E 
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Figure B.1l. Spherical Trigonometry Reference 


cos 8 = - cosa - cos y + Sina- sin y + cos B/E 


where E is the Farth's radius. 

If in this coordinate system a target has position 
(latitude, longitude,altitude) = (X,Y, ,A,) and Great 
emrcle course Cy at time 1, and if by time 2 it has moved D 
units at constant altitude, the new position (X,,Y5,A,) 


and course Cc, are computed by: 


7 ei . 
Yo =—esin: » (Sin Yy - cos D+ cos Y,- sin D -cos C,) 


X — oe she teke D = Sin C€ 


D l ~ cos Yo) 


1 


= _ -l,_. _ a ae 7 : _ 
Cc =——eCt— Gos —-(Ssin Cy sin(X, X,) Sime vy cos Cy cos (X, X,) 


x (C) + 1C,)) 


i 








XK, 0) (x,, 0) 


Figure B.2. Spherical Movement 


and 
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ABP SND EX Cc 


FIGHTER MANEUVERS 


Only one set of fighter maneuvers are analyzed in this 
study. Some types of maneuvers are certainly better than 
others and the EKF algorithm has been demonstrated to be 
robust enough to operate well with a wide variety of 
different maneuvers; the essential features are that there 
be some variation in motion parameters [Ref. 9] and that the 
fighter be able to determine with reasonable precision his 
Own position. 

The fighter maneuvers used for these simulations are as 
follows: 

(a) Upon initial target observation he immediately 
(instantaneously) sets a course 0.7 radians (about 40 
degrees) left of the first bearing observation. 

(b) The fighter maintains a constant course (constant 
speed and altitude with no horizontal turns--in the spheri- 
cal model this is a Great Circle track) until either a 
pre-set number of new observations have been taken, or until 
the relative bearing angle to the target exceeds 0.96 
radians (about 55 degrees) right. At this time he turns 
(again instantaneously) to establish a new course 0.7 
radians to the right of the most recent bearing observation, 


and similarly continues until relative bearing exceeds 0.96 
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radians left or until a set number of observations have been 
taken. 

(c) This not-quite-symmetrical Zig-zag pattern continues 
until a set number of total observations have been recorded 
or until the target has passed the fighter (bearing obser- 


vation angle less than - > or greater than =). 
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APPENDIX D 


ELEVATION ANGLE CORRECTION FOR EARTH'S CURVATURE 


The Extended Kalman Filter algorithms proposed here were 


developed based on a "flat-Earth" motion model. To enhance 


the performance of the bearing-elevation EKF when using 


inputs from a spherical motion model, a simple, approximate 


elevation angle correction was developed. 


Referring to the figure on the next page, the notation 


is as follows: 


oe; = 


Then 


Ir 


-a) = sin(>-8) = cos g = [ 


true elevation angle measured from the 
Heri Zon 


"corrected" elevation angle to be input to 
the EKF 


line-of-sight range 


radius of the Earth (assumed here to be 
3444nm) 


altitude of target 


altitude of fighter 


F 
=~ + 
a 


At 


~ Ap 
ee by Construction 
ie og , | 
a aa by small angle sine approximation 
2 2 2 
R +(E+A,) ~ (E+A,) 
Zo eA) 


[ 


] 
F 


by law of cosines 


2 2 Z 
R +(E+A,) - (E+A,,) 


2 oR (E+A,) 


] by small angle sine 


approximation 
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Figure D.1i. Elevation Angle Correction 


A -Ap R°+(E+A_) “-(E+tA_) * 
eo ee 
R 2 R (E + A,) 





| 


‘ 
3 
Ile 


R°-AC+2A,A “AS 
2 R (E+A,) 


SO, since A. and A, are "small" and their difference is 


somal 1" 


ra a 
Ze 
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t= a + Be is the correction (based on EKF-estimated 


SO @ oe 


range). 
NOTE: As the magnitude of An,-Ap increases relative to 


range, the approximation becomes less accurate, but as 


range decreases the correction goes to zero. 
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APPENDIX E 


APL COMPUTER PROGRAM 


The APL program shown here assumes global variables 
as defined in the initial comment section. The program 
updates fighter and target positions as described in 
Appendices B and C, generates random observations, and 
computes corresponding BO and BE EKF estimates. Choice of 
output variables and programming of output statistics is 


left to the user. 
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